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DEVELOPABLE  WINGS* 


/"Following  is  a  translation  of  an  article 
written  by  S.  V.  Vellander  in  Vestnlk  Le- 
nlngradskoRO  Unlverslteta  (Herald  of  Lenin¬ 
grad  University),  No. 197  Ser.  Mat.  Mekb.  i 
Astron.,  No.  4,  1959  >  pages  113-120._7 

In  this  work,  we  shall  develop  a  method  of  computing 
aerodynamic  characteristics  of  wings  of  finite  span,  of  a 
type  we  shall  call  developable. 

The  whole  study  shall  proceed  under  the  assumption 
that  the  gas  is  devoid  of  forces  of  internal  friction  and 
moves  adiabatic ally.  We  shall  assume  that  the  strong  dis¬ 
ruption  surfaces  do  not  recoil  from  either  the  frontal  or 
the  rear  edges  of -the  wing.  The  motion  shall  be  assumed 
to  be  fixed.  The  interference  of  the  wing  with  the  fusela¬ 
ge  shall  be  left  out  of  the  examination. 


1 .  Developable  Wings 


We  shall  consider  wing  finite  span  of  EFGH,  stream* 
lined  by  the  supersonic  flow 
of  the  gas,  to  be  developable 
if: 

1 )  the  upper  and  lower  sur¬ 
faces  of  the  wing  are  develo¬ 
pable  surfaces; 

2)  the  frontal  edge  of  the 
wing  (EH)  is  a  straight  line; 

3)  the  rear  edge  of  the  wing 
(FG)  and  the  consolic  edge  (EF) 
are  such  that  any  Mach  line- 
traced  forward  from  their  points 
intersect  the  airfoil  only  at 
points  on  the  frontal  edge. 


2.  Generalization  of  the  Prandl-Mayer  Flow 
xPlease  see  page  -a- 
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The  following  feet  is  well  known  from  geometry?  for 
any  developable  surfaces  (except  cones  and  cylinders) ,  a 
curve  of  -  double  curvature  can  be  found  having  the  property 
that  the  envelope  of  the  tangent  planes  to  the  curve  will 
yield  the  initial  developable  surface  S. 

This  purely  geometrical  fact  allows  for  an  easy  ge¬ 
neralization’  of  the  "flows  of  Prandtl-Mayer,  and  enables  us 
to . obtain  analogous  flows  about  any  arbitrarily  chosen, 
convex,  developable  surface  —  including  cones  and  cylin¬ 
ders. 

Indeed,  let  us  examine 
an  arbitrary  developable  sur¬ 
face  (Fig.  2).  We  find  on  the 
basis  of  the  geometrical  fact 
cited  that  it  can  be  approxi¬ 
mated  to  within  any  degree  of 
exactness  by  a  surface  compo¬ 
sed  of  flat  strips,  provided 
the  number  of  such  strips  la 
sufficiently  large.  Cones  and 
cylinders  shall  be  temporarily 
excluded  from  consideration.  The 
possibility  of  such  an  approxi¬ 
mation  for  cones  and  cylinders  does  not  follow  from  the  ge¬ 
ometrical  fact  stated  but  is  directly  evident. 

Thus,  we  may  maintain  that  any  developable  surface 
can  be  approximated  by  a  surface  composed  of  flat  strips, 
provided  their  number  is  sufficiently  large. 

We  shall  consider  a  convex  developable  surface  which 
is  flat  to  the  right  from  the  rectilinear  generatrix  (AB) 
and  curvilinear  to  the  left  of  same.  We  shall  assume  that 
a  gas  flows  along  the  flat  part  of  this  surface  with  super¬ 
sonic  velocity  V]  with  constant  hydrodynamic  elements.  The 
supersonic  velocity  should  be  such  that  the  projection 
of  Vj  on  the  normal  to  the  generatrix  AB  traced  in  the  pia¬ 
no  ABCB  remains  supersonic. 

We  shall  try  to  solve  the  problem  of  the  movement  of 
the  gas  about  the  curvilinear  part  of  the  developable  sur¬ 
face  under-  consideration. 

For  this  purpose,  we  shall  approximate  the  develo¬ 
pable  surface  by  a  surface  consisting  of  flat  strips.  The 
process  of  solving  the  problem  becomes  quite  clear. 

Indeed,  the  streamlining  mechanism  of  a  surface  com¬ 
posed  of  flat  strips  (Fig,  2)  is  very  simple.  The  hydrody¬ 
namic  elements  on  such  strips  are  constant,  while  at  the  ed¬ 
ges  AB,  A |  B?  i  AgB2. . . ,  we  always  have,  the  well-known  sli¬ 
ding  streamlining  over  an  evolving  angle.  Here  the  projec¬ 
tion  of  the  velocity  upon  the  edge  of  the  angle  remains- 
constant,  and  has  no  meaning  for  streamlining. 
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Let  us  trace  a  plane  It  along  such  an  edge  and  the 

velocity  vector,  , 

We  denote  'by  Kc  the  projection  of  the  velocity  on  the 
d  r.f  the  @d«e.  and  by  t4i  the  projection  of  the  velo-* 

c it^ ^ on°t he  normal  to  the  edge  that  lies  in  the  plane  II. 

Me  shall  further  provide  with  the  index  1  the  va¬ 
lues  of  all  hydrodynamic  elements  prior  to  streamlining  &- 
cross  the  angle,  and  those  streamlining  past  with  the  in¬ 
dex  g.  The  angle  of  rotation  of  the  gas  flow  shall  be  de¬ 
noted'  bv  A  f%  .  Then,  to  within  errors  amounting  only  to  the 
of  small  magnitudes,  we  shall  have 

V\e5 

¥*,  a  VV  t  ~h  A  p,  (2,  S) 

y v*.  t  ~  a 


higher  power 


Hvl+vS  +  i—i — 


where  «&•  is  the  velocity  of  sound,  t#  the  constant  in  Eer- 
boulll’s  formula  and  Jt  the  exponent  of  the  ad ia antics. 

A  subsequent  application  of  the  equations  (2,1)  and 
a  corresponding  sequence  of  projections  on  the  changing  ax- 
rections  of  the  axes  H  and  permit  to  compute  the  g&s 
flow  to  within  any  degree  of  precision  desired. 

However1  s  it  is  mono  convenient  to  oh^&.xn  the  onuixin*" 

r-y  differential  equations  for  \/x  and  <4i*  „ 

*  To  derive  these  e- 

quations  we  shall  examine 
Fig,  3.  Let  the  edge  AB 
be  subject  to  a  flow  pos¬ 
sessing  velocities  riij  and 
/xj  ,  Let  us  suppose  that 
the  angle  of  rotation  of 
the  flow  about  the  edge  AB 
is  Ah  and  let  the  angle 
between  the  edges  AB  and 
A‘B?  be  A^.  Then  for  the 
velocities  Vh,3  and  Y  w,  3 
with  which  the  gas  flow  ap. 

proaches  the  edge  A*B*  we  shall  have  the  equations* 


^*3*  Vn,iCQ&  Ay  —  l/r,j.slnAy, 
Vt,3~  V^cosAy-f  y^xSinAj. 


M 


On  substituting  V^A  and  V<,  &  from  (2,1}  .into  (2»2)  ®.nd  re 
placing  sines  by  angles  and  cosines  by  units,  we  oowxn, 

v. 

Vvt  ,-rf 


n,  3 


t.  I 


Ay. 


Vtj  j —  V'*  i  — 


(2,3} 
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On  dividing  poth  parts  of  the  equalities  (2,3)  by  &n& 
letting -a/3  approach  zero,  we  obtain  the  differential  equa¬ 
tions  required. 

‘iX* *X&  v  az 
**  '  \fv%~*.  u as  * 

d/r  _  v  d~*  (2,4) 


To  these  equations  we  must  still  adjoin  Bernoulli’s 
integrals 

V£~t~  V?)  ~j~  /^~Ti  x  h- 

If  the  equation  of  the  developable  surface  is  given, 
then  the  function  is  a  completely  defined  function 

of  the  variable/?  *  Moreover,  the  parameters  of  the  gas 
flowing  along  the  flat  part  of  the  surface  provide  us  with 
the  initial  conditions.  Hence,  the  system  (2,4)  can  be  in¬ 
tegrated  on  taking  Into  account  (2,5),  and  will  determine 
the  completely  defined  functions  V-n  and  %  ,  which  depend 
on  the  argument/3.  After  the  Integration  of  (2,4),  on  ta¬ 
king  into  account  (2,5)  with  the  known  V-ti  and  Mrr,  we  shall 
f  ind  f  and/5,  adjoining  to  the  equation  (2,5)  the  adiabatic 
condition  p 

-/T  *  const.  (2,6) 

This  solves  the  above  problem. 

3.  Application  of  the  Generalized  Flows  to  Developable 
Wings 

Let  EFGH  be  the  developable  surface  in  Fig.  1 ,  and 
EH  shall  be  the  wing’s  frontal  edge.  Let  us  suppose  that 
the  wing  meets  the  oncoming  gas  with  the  velocity  Vj*  De¬ 
pending  on  the  relative  position  of  the  wing  with  respect 
to  the  oncoming  gas,  one  or  two  surfaces  of  strong  disrup¬ 
tion  are  formed  about  the  edge  EH  (above  and  below),  but 
only  one  sliding  streamlined  flow  at  the  angle  of  the  ro¬ 
tation. 

Generally,  the  surfaces  of  strong  disruption  will 
not  be  planes.  However,  with  properly  chosen  (and  thus  the 
only  ones  of  Interest)  profiles  of  the  wings  the  curvatures 
of  these  surfaces  can  be  neglected;  the  resulting  precision 
corresponding  to  the  third  approximation  of  the  theory  of 
the  thin  wing.  We  can  also  assume  that  beyond  the  surfaces 
of  strong  disruption  the  entropy  is  constant  and  whirls  are 
absent. 

The  argumentation  of  these  statements  does  not  differ 
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in  principle  from  that  used  in  the  gas  dynamics  of  planar 

flows:  the  computation  process  is  only  slightly  more  com- 
plicated*  With  whirls  absent  and  the  entropy  constant  (be¬ 
yond  the  surfaces  of  strong  disruption)  we  shall  have  the 
generalized  flow  of  Prandtl-Mayer  with  the  cited  precision. 

The  above  Is  quite  sufficient  for  evolving  a  scheme 
for  the  determination  of  the  aerodynamic  characteristics 
of  the  developable  wing. 

First,  we  must  consider  the  sliding  streamlining  pro¬ 
blem  of  the  frontal  edge.  From  the  solution  we  shall  get 
4  4.8  and  ft.  when  the  flow  is  parallel  to  the  surface  of 
the  wing  at  the*  frontal  edge.  Then,  given  the  Initial  con¬ 
ditions  Kr*  *  etc we  must  integrate  the  system  (2,4)  and 
solve  for  vn~  *«<?).  V,  -  V,  fp).  p**p(})  r  p  =  fMp). 

And,  while  working  on  the  solution  of  the  problem  of  the 
sliding  flow  over  the  frontal  edge  of  the  wing,  one  must 
keep  in  mind  that  remains  constant  both  during  the  tran¬ 
sition  over  the  disruption  and,  in  that  Prandtl-Mayer  flow 
which  may  result  at  the  frontal  edge, 

4,  Certain  Remarks 


Note  1.  In  a  large  number  of  cases,  the  integra¬ 
tion  of  the  system  (2,4)  can  be  performed  in  series.  One 
has  only  to  express  and  V+  as  series  in  /?  and  to 

obtain  from  system  (2,4)  the  coefficients  of  the  series  for 
/■»  and  • 

TXWW  is  expressed  in  the  form  of  a  series,  same 
results  can  be  obtained  by  the  method  of  successive  appro¬ 
ximations. 

Should  we  use  a  series  expansion  for  the  flow  about 
the  wing  in  the  presence  of  surfaces  of  strong  disruption, 
we  must  delete  from  that  series  terms  with  higher  powers  of 
(\  than  the  third;  the  errors  will  be  of  order^/?^ 


Note  2,  Should,  the  streamlined  developing  surface 
be  the  surface  of  a  circular  cone,  the  quantity  Jjr/d/i 
will  be  constant  and  the  integration  of  the  system  (2,4)  Is 
simplified. 

Indeed,  let  &  _  ,, 

S ”~c.  (4,1) 


Then,  dividing  the  first  equation  of  (2,4) 
we  have  the  single  differential  equation; 

wv*  i  &  • 

C  V* 

Integrating  by  quadratures  we 

from  the  second  equation  of ,l 2,4), 

♦  V«—  Vm.9,  V4*  Vt,*  />«*/>,  m  p«“*p0, 


by  the  second. 


(4,2) 

shall  get 
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Note  3,  The  formulas  become  especially  simple  when 

the  streamlined  surface  is  a  circular  cone,  and  If  the  pro¬ 
file  is  so  thin  that  first  powers  of  fl  are  sufficient*  In 
this  case,  all  computations  become  easy. 


Note  4.  In  the  case  of  the  circular  cone  the  magni¬ 
tude  Jjf/JfS  has  a  very  simple  geometrical  meaning. 


If  M  is  a  point  on  the  surface  of  the  cone,  L  its 
distance  from  the  apex,  and  R  the  radius  of  curvature  at  M 
for  the  curve  formed,  toy  the  intersection  of  the  cone  with 
a  plane  passing  through  M  and  normal  to  the  generatrix,  then 


5.  Flow  Equations  for  Conic  Wings 


Of  all  developable  wings  the  conic  ones  are  of  the 
greatest  interest*  It  is  therefore  greatly  desirable,  at 
least  for  such  wings,  to  obtain  a  transformation  of  the  e- 
quations  in  (2, 4)  such  as  would  explicitly  introduce  into 
these  'equations  the  characteristics  of  the  conic  wing's  pro¬ 
file,  This  transformation  will  be  treated  in  this  fifth 
section  of  the  author's  work. 


Let  us  suppose  that  A BCD 
is  the  profile  of  a  conic  wing 
of  the  points  of  contact  with 
the  fuselage  (Pig,  4) ,  We 
shall  choose  for  the  coordina¬ 
te  plane  the  vertical  pla¬ 
ne  containing  the  profile  ABCD* 
We  shall  direct  the  x  axis  in 
the  direction  of  the  unpertur¬ 
bed  flow,  the  y  axis  vertical¬ 
ly  upward,  and  the  s  axis  per¬ 
pendicular  to  x  and  y. 

Let  us  suppose  that  in  the 
chosen  coordinate  system  the  a- 
pex  of  the  cone  E  has  the  coor¬ 
dinates: 


*  **  Q.  y  ™  O,  z—  l 


0J) 


whether  this  apex  does  or  does  not  actually  belong  to  the 
wing. 

And  let  us  suppose  further  that  the  profile  ABCD  is 
given  by  equations 

2*0  (5,2) 
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in  the  chosen  coordinate  system. 

Then,  for  every  value  of  the  parameter  t  shall  cor¬ 
respond  a  generatrix  L  of  the  Cone,  given  by  the  equations: 

Jt  9  X~~  l  /c?\ 

~T  •  (§3> 

The  directional  cosines  of  this  generatrix  L  shall 

be  obviously  given  by  formulas:  v 

COS {L,  X)  S*  —  -rr«s=======:  , 


A 

cos  (L,  y)  «* 
cos  (L?z)  ■ 


l 

yw+i^Tb- ' 


In  addition  to  the  straight  line  L  we  shall  examine 
the  tangent  T  to  the  airfoil- ABOD.  This  line  is  given  by 
the  equations: 

y-?(0-V(0<*-0. 

*  *=•  0  <5,3> 


and  It  has  the  following  directional  cosines: 

COS  (  7*fx)  —  -  ■■jsasaaaa.  , 

V  Vl  +  T'* 

COS  (  Ky)  «  ■  , 

V  1  +  f  * 

cos  (T,  z)~  0. 

From  the  formula: 

cos (f^L)^  cos {TTx) cos  (L.  x)±cos(T,y)coi{L,  y)  f- 

4-  COS(r.  2)  COS  (I,  2) 

and  from  formulas  (5,4)  and  (5,6)  we  obtain  directly: 

,  co*(T'L)~-7r{7^^ 

and  consequently :  r  +»  y 


sin  ( r,  £.) 


VW-9)*  +  <M»  +>*> 

V'l  4  ?  *  f1  +  ?*  +  ?* 
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♦ 


« 


We  shall  examine  on  the  profile  ABCD  two  points  cor¬ 
responding  to  the  parametric  values  t  and  t  +  dt. 

The  angle  between  the  generatrices  passing  through 
these  points  shall  be  the  angle  between  the  normals  to 
the  surface  of  the  wings  for  these  points  shall  be  #*/?• 

We  shall  therefore  determine  the  two  magnitudes 
and  dy  ***** 

If  we  denote  the  distance  of  a  given  point  on  the 
profile  from  the  apex  E  of  the  cone  by  R  and  the  distance 
between  near  points  on  the  same  profile  by  dS,  we  Shall 
ha we  %  A 

d^^JJxJ:lds.  (5,U>) 

Since  the  formula  for  the  sine  (T?k)  is  known,  and 

(M» 


(5,12) 


(5,13) 


since  R  and  dS  are  given  by  formulas! 


we  obtain 


R  -  Vi-  +  f*  +  t\ 

ds  ~  YTTVldt, 

dyt  **  (0  dtf 


where  we  let 


.  ...  W -  oi’-f  />'0  ±  v')'1 

- *8..*  4  it. 


We  shall  now  determine  the  second  magnitude  of  inte¬ 
rest,  namely  d/». 

We  introduce  the  unit  vectors 

Lt **  icu%{L,  x)  -f*j  £q£>(_L,  y)  ~t  k,co&{L,  z% 

f,  ~  Tcos  ( 7pr)  -tjcos  ( rfy)  •+  k  COS  ( (5, 1 4) 

where  1,  j  and  1c  are  the  orthonormal  of  the  coordinate  axes. 
It  is  obvious  that  the  vector 


N 


1  =  Ti  X  Lt 


(5,15) 


has  the  direction  of  the  inner  normal  to  the  surface  of  the 
wing. 

By  using  the  known  formulas  for  vector  products  and. 
the  previously  found  expressions  for  the  direction  cosines, 

we  shall  have  . . — « — - 

*•+  ixR  )  (5  16) 


N  sin  ( r,  L)  ** 


Ny 


(i 

y 


vT+*a 


(5,12) 
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FT??*' 

f?'  —  y  < 

K"i  -+-  #'*■ 


(5,18) 

($19) 


« 


Therefore  ^ 
cos  (A/,  x) 


4''*V 

cos (iV,  y) 


vV*'™?>M-/A(i~Hrs) 

_ _ ~l 

y  (ty'— +  <?'■*} 


/*<*>, 

A(*), 


(5,20) 


COS  ~ 


FeF 


-f* 


/»<<), 


where  the  functions  f t ,  fg»  and  f,  are  merely  abbreviated 
notations  for  the  magnitudes  in  the  left  members* 

If  the  formulas  of  (5,20)  determine  the  direction 
cosines  of  the  normal  corresponding  to  the  parameter  L,  it 
is  obvious  that  the  magnitudes 


Ax  =/,  in  +f>  {Qdt, 

«  Ait)  +■  (5,20 

Az*~Mt)  +  fkt)dt 


will  be  the  directional  cosines  of  the  normal  corresponding 
to  the  value  of  the  parameter  t  +  dt. 

On  considering  the  unit  vectors 

jlcev  ■««,  v,r*iL  <ei*' 

t  ®  if* +  #-  +  *&.  (s,22l 

A  » iAf  -yjAy  ~t  kAjr, 


we  find  easily  that 

d§  «*  i  » j  j  A  i  sin  #  =»  i«X  A{. 

with  a  precision  to  the  small  values  of  the  higher  orders. 

On  composing  the  modulus  of  the  vector  product  of 
n  and  A,  we  obtain 


in  which  we  let 


*>  4>a  (*)<#, 


(5,23) 


'.*.5. 


US*. 


V  Ifjf*  ~  fufit  +  (A /» —  /i/af  +  {/ A *~A/<)  ■  £5,24} 


Formulas  (5,12)  and  (5,23)  enable  us  to  pass  in  the 
equations  (2,4)  from  the  independent  variable  Z3  to  the  inde¬ 
pendent  variable  t. 

On  replacing  4/S  and 4jr in  these  equations  by  their  ex¬ 
pressions  according  to  formulas  (5,12)  and  (5,23) »  we 
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■> 


\ 


the  equations 


r~%^  -  w«t>»  (#), 

it  V  v% 


iV-c 

"37 


vvm% 


which  solve  the  streamlining  problem  of  the  wing,  whose 
basic  principle  is  given  by  the  equations  (5,2), 

On  integrating  the  equations  (5, 25) ,  we  find  Vf?  and 
as  functions  of  the  argument  t. 

After  that,  the  pressure  p  also  will  be  a  known 
function  of  argument  t. 

With  p  known,  and  with  expressions  available  for  the 
directional  cosines  of  the  normal  as  well  as  for  the  magni¬ 


tude  //;  the  final  computation  of  the  aerodynamic  factors 
presents  no  difficulties. 

An  analogous  transformation  can  be  made  for  any  ar¬ 
bitrary,  developable  surface. 


v 
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END 


2023 


10 


